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ABSTRACT  :  We  are  investigating  the  fine-scale  structure  of  the  Earth's  crust  in  order  to 
understand  its  influence  on  high  frequency  (0.5-30  Hz)  regional  wave  propagation.  We  have 
developed  geologically  based  stochastic  models  of  crustal  heterogeneity  for  a  number  of  different 
tectonic  provinces  and  crustal  levels.  Here  we  have  focused  on  examining  geologic  maps  from 
extended  terranes,  and  seismic  data  from  the  northern  and  southern  Basin  and  Range.  We  have 
developed  two  different  deterministic/stochastic  velocity  models  which  produce  synthetic 
seismograms  showing  many  of  the  reflectivity  patterns  and  crustal  events  (Pg,  PmP,  Sg,  SmS) 
observed  in  Basin  and  Range  near-vertical  and  wide-angle  seismic  data.  We  also  have  developed 
several  ways  to  measure  the  lateral  heterogeneity  scale  length  from  unstacked  seismic  data,  and 
have  tested  them  on  synthetic  data  with  promising  results. 
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OBJECTIVE  :  We  are  investigating  the  fine-scale  seismic  velocity  heterogeneity  of  the  Earth's 
crust  over  a  broad  range  of  scales  (-0.025  to  ~35kin)  in  an  effort  to  understand  how  wavelength 
scale  heterogeneity  affects  seismic  wave  propagation  in  the  the  cmst  in  the  0.5-30Hz  band.  We  are 
approaching  this  problem  in  three  ways.  1)  Estimates  of  crustal  velocity  heterogeneity  can  be  made 
from  geologic  maps  and  petrophysic^  data,  and  from  seismic  reflection  data.  Statistical  analysis  of 
geologic  maps  and  petrophysical  data,  from  which  stochastic  models  for  velocity  heterogeneity  are 
constructed.  The  stochastic  models  of  different  types  of  lithologies  are  then  used  to  develop  whole 
crustal  models  for  which  synthetic  seismograms  can  be  calculated.  To  improve  the  direct  image  of 
crustal  velocity  heterogeneity  in  the  crustal  column  we  are  reprocessing  COCORP  reflection  data 
from  the  Basin  and  Range.  2)  Synthetic  seismograms  calculated  for  the  stochastic  models  are 
compared  to  crustal  reflection  and  refraction  seismograms.  3)  We  are  evaluating  methods  for 
estimating  lateral  heterogeneity  scale  lengths  using  several  correlation  measures  on  unstacked 
synthetic  seismic  data. 

RESEARCH  ACCOMPLISHED  :  Reprocessing  the  COCORP  4(PN  Transect  -  Nevada  Line  2. 
Vertical  incidence  seismic  data  from  the  northern  Basin  and  Range  Province  are  often  characterized 
by  a  transpiirent  upper  crust,  a  highly  reflective  midcrust,  a  slightly  less  reflective  lower  crust,  and 
a  strong  reflection  Moho  up  to  a  second  in  length.  The  COCORP  40ON  transect  -  Nevada  line  2  is 
part  of  a  vibrator  seismic  survey  that  extends  across  the  western  Cordillera  of  the  United  States. 
Line  2,  situated  in  the  Basin  and  Range  Province,  has  been  processed  and  interpreted  previously 
(Klemperer  and  others,  1986;  Allmendinger  and  others,  1987:  Hauge  and  others,  1987).  A  joint 
vertical-incidence  to  wide-angle  reflection/refraction  explosion  experiment  (known  as  the 
PASSCAL  Basin  and  Range  Experiment)  was  conducted  over  the  same  area  in  1986. 

Interpretations  of  these  two  datasets  have  produced  somewhat  contradictory  models  of  the 
cmst  in  the  Basin  and  Range.  One  COCORP  CMP  section  is  a  single  fold  near-trace  CMP  section, 
which  shows  a  pair  of  bright  reflectors  at  the  base  of  the  cmst,  whereas  the  stacked  CMP  stacked 
section  shows  a  broad  zone  of  reflectors  from  9.5  -  10.5  s  two-way-time  (TWT),  as  opposed  to 
two  distinct  bands  of  reflectivity.  Neither  section  shows  bright  or  continuous  middle  crustal 
reflectivity  across  the  section,  whereas  a  limited  (~20km)  dataset  of  densely  recorded  near  vertical 
incidence  explosion  data  shows  pronounced  reflectivity  in  the  middle  cmst  (at  4-7  s  twt).  The 
disparity  between  the  stacked  and  single-fold  vibrator  and  explosion  images  of  the  middle  to  lower 
cmstal  reflectors  may  be  caused  by  inadequate  seismic  processing  of  the  COCORP  data.  Surface 
related  problems  were  pointed  out  by  the  original  processors  (Klemperer  et  al.,  1986).  It  is  our 
intention  to  use  current  processing  capabilities  to  remove  the  near-surface  effects  and  produce  a 
better  stacked  image.  Getting  reliable  consistent  images  of  the  middle  cmst,  the  Moho  transition 
zone  is  cmcial  to  estimating  the  parameters  needed  for  stochastic  models.  We  will  then  have  a 
reflectivity/velocity  model  for  the  cmst  and  mantle  in  an  extended  terrane  which  can  be  used  to 
generate  synthetic  seismograms. 

To  reprocess  COCORP  Nevada  Line  2  we  followed  the  processing  flow  for  Nevada  line  2 
outlined  in  Allmendinger  and  others,  1987,  to  reproduce  the  original  COCORP  image  and  give  us 
a  foundation  on  which  to  make  improvements.  Additional  trace  editing  and  careful  muting 
removed  some  of  the  ground  roll  and  extraneous  noise  later  in  the  section.  Velocity  analysis  was 
done  interactively  with  a  starting  velocity  model  taken  from  the  1986  PASSCAL  wide-angle  survey 
of  the  same  region  (Holbrook,  1990).  The  stack  produced  after  the  velocity  analysis  showed 
incremental  improvements.  We  further  improved  the  section  by  calculating  and  applying 
refraction  statics,  which  strip  away  some  of  the  effects  of  the  irregular  topography  and  surface 
layers  in  the  basins  and  ranges.  The  CMP  stack  following  refraction  statics  and  several  iterations 
of  velocity  analysis  show  continuous  Moho  and  middle  crustal  reflectors  across  the  section. 
Additional  processing  will  include  residual  statics  and  time  and  in  depth  migration. 

The  reprocessed  image  (Figure  1)  appears  much  more  consistent  with  the  patterns  seen  in 
the  limited  explosion  data,  and  thus  suggests  that  we  may  be  able  to  process  the  entire  COCORP 
Basin  and  Range  dataset  to  estimate  crustal  fabric  across  a  regional  scale. 
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Construction  of  stochastic  model  for  the  northern  Basin  and  Range.  Figure  2  displays  our  velocity 
model  for  the  northern  Basin  and  Range.  The  transparent  upper  crust  is  modeled  as  felsic 
intrusions  into  a  gneissic  basement  terrane.  The  horizont^  and  vertical  ch2iracteristic  scales  for  the 
upper  crustal  zone  are  963  m  and  817  m  respectively,  with  a  fractal  dimension  of  2.5,  and  a  binary 
probability  density  function  (pdf)  of  67%  fast  (avg  6.0  km/s)  and  33%  slow  (avg  5.7  km/s). 
These  statistics  were  obtained  from  outcrop  maps  within  the  Chocolate  Mountains  of  southeastern 
California.  (Although  not  within  the  northern  Basin  and  Range  Province,  maps  of  this  area  show 
the  effects  of  upper  crustal  Mesozoic  plutonism).  Seismic  velocities  are  taken  from  laboratory 
measurements  of  rocks  with  similar  compositions  and  are  consistent  with  the  bulk  velocities 
determined  from  upper  crustal  refractions. 

High  heat  flow  and  historical  volcanism  throughout  the  Basin  and  Range  indicate  that  mafic 
intrusions  are  likely  to  exist  at  several  levels  in  the  crustal  colunm.  These  intrusions  provide  the 
impedance  contrasts  required  to  produce  high  amplitude  reflections.  Rheological  arguments 
suggest  that  these  intrusions  should  have  different  orientations  depending  on  their  depth  of 
emplacement,  with  horizontal  intrusions  prefering  weak  rheological  boundaries,  such  as  at  the 
brittle/ductile  transition  and  at  the  crust/mantle  boundary  (Parsons  et  al.,  1992;  Holliger  and 
Levander,  1994).  For  these  reasons  we  have  developed  a  stochastic  model  of  velocity 
heterogeneity  consisting  of  either  horizontally  or  vertically  elongated  high  velocity  bodies  within  a 
more  felsic  matrix.  The  zones  having  elongated  horizontal  bodies  within  the  midcrust  and  at  the 
Moho  produce  the  observed  high  amplitude  reflectivity.  The  mid  and  lower  crust  consists  of  dikes 
and  sills  with  characteristic  scales  of  200  m  and  800  m  depending  on  orientation,  a  fractal 
dimension  of  2.7,  and  a  binary  pdf  of  5%  fast  (avg  6.8  km/s)  and  95  %  slow  (avg  6.5  km/s).  The 
fractal  dimension  is  consistent  with  that  measured  from  outcrop  maps  of  the  Ivrea  Zone  in  northern 
Italy,  an  area  believed  to  expose  extended  and  intruded  lower  crust  similar  to  that  in  the  Basin  and 
Range.  Though  we  presently  have  no  seismic  constraints  on  the  percentage  of  intrusion  in  the  pdf, 
we  have  chosen  a  small  extreme  to  emphasize  the  importance  of  orientation  in  producing 
reflectivity.  Average  crustal  velocities  and  gradients  were  obtained  by  averaging  the  velocity 
models  shown  in  Figure  3. 

Construction  of  stochastic  model  for  the  southern  Basin  and  Range  in  southeastern  California.  In 
contrast  to  the  highly  reflective  crust  prevalent  in  the  northern  Basin  and  Range,  southeastern 
California  north  of  the  Chocolate  Mountains  exhibits  a  comparatively  nonreflective  seismic 
signature.  We  have  therefore  modeled  this  area  without  mafic  intrusions,  and  instead  have  used 
the  statistics  derived  from  maps  of  Mesozoic  plutons  intruded  into  Proterozoic  North  American 
crust  in  the  Chocolate  Mountains  (Figure  2,  upper  crust)  as  describing  the  entire  crustal  column. 
Average  crustal  velocities  and  gradients  were  obtained  from  the  analysis  of  the  PACE  92  refraction 
profile  of  this  same  area  (Parsons  and  McCarthy,  submitted  to  Journal  of  Tectonics). 

Modeling  PASSCAL  86  northern  Basin  and  Range  data  and  RISC  92  southern  Basin  and  Range 
data.  The  two  velocity  models  described  above  were  used  to  generate  2-D  visco-elastic  finite- 
difference  synthetic  seismic  sections.  For  the  vertical  incidence  simulations,  the  models  were  20 
by  36km  with  20  m  grid  spacing.  The  source  pulse  was  a  Ricker  wavelet  with  an  18  Hz  center 
frequency.  For  the  wide-angle  simulations,  the  models  were  210  by  36  km  with  60  m  grid 
spacing.  The  source  was  a  Ricker  wavelet  with  a  5  Hz  center  frequency.  The  near-surface  was 
modeled  as  a  series  of  2-D  layered  sedimentary  basins.  Average  velocity  within  the  near  surface 
layer  was  3. 1  km/s  with  a  gradient  of  0.8  km/s/km  and  a  Q  of  50.  Velocities  were  determined  from 
travel  time  analysis;  Q  was  determined  from  spectral  ratios  of  the  near-offset  first  arrival. 
Sedimentary  layers  were  modeled  as  fractal  tabular  bodies  with  horizontal  and  vertical  characteristic 
scales  of  6000  m  and  60  m,  respectively,  a  fractal  dimension  of  2.5  (the  smoothest  possible  binary 
fractal  fabric),  and  a  binary  pdf  of  50-50%  high  and  low  velocities  with  velocity  contrasts  of  0.4 
km/s. 

The  results  from  the  PASSCAL  86  northern  Basin  and  Range  simulations  are  shown  in 
Figure  4.  Figure  4a  shows  a  shot  record  (Shotpoint  4B)  from  the  1986  PASSCAL  seismic 
experiment,  and  Figure  4b  shows  our  vertical  incidence  synthetic.  Both  show  a  reflective  midcrust 
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and  a  1  second  thick  Moho.  No  bulk  changes  in  velocity  occur  at  the  onset  of  the  reflective 
midcrust. 

The  results  from  the  RISC  92  southern  Basin  and  Range  simulations  are  shown  in  Figures 
5.  The  seismic  section  displayed  in  Figure  5a  is  a  composite  gather  of  a  number  of  shots  into  a 
stationary  recording  spread  near  the  Chocolate  Mountains  of  southeastern  California.  The 
reverberatory  nature  of  wide-angle  PmP  is  reproduced  remarkably  well  in  the  synthetics  (Figure 
5b)  by  a  combination  of  crustal  scattering  and  basin  reverberation.  Coda  amplitude  analysis  has 
shown  that  the  initial  coda  is  dominated  by  near  surface  reverberation,  but  the  later  coda  is  due  to 
body  wave  scattering  and  phase  conversions.  Crustal  heterogeneity  is  also  necessary  to  produce 
the  strong  SmS  arrival  observed  in  the  data. 


Statistical  measures  of  crustal  heterogeneity  :  The  Basin  and  Range  data  show  that  describing 
highly  reflective  zones  within  the  crust  as  stochastic  media  can  produce  synthetic  seismograms 
having  many  of  the  features  seen  in  seismic  reflection  and  refraction  data.  These  crustal  models  are 
now  based  on  the  statistics  derived  from  digitized  maps  of  middle  and  lower  crustal  exposures. 
Synthetic  seismograms  generated  through  these  models  show  that  their  fine-scaled  features  affect 
the  propagation  of  seismic  signals  at  and  above  0.5  Hz  frequency  (Holliger  and  Levander,  1992; 
Holliger  et  al.,  1993;  Levander  et  al.,  1993). 

Our  objective  is  to  relate  the  statistics  of  the  scattered  wavefield,  observed  in  the  synthetic 
seismograms,  to  the  statistics  of  the  underlying  medium.  In  the  process,  we  explore  the  parameter 
space  that  determines  the  statistics  of  the  medium.  This  includes  modality,  fractal  dimension, 
velocity  fluctuation,  fabric  anisotropy,  and  structure  of  the  velocity  pdf.  The  current  work  focuses 
on  relating  the  statistics  of  the  wavefield  to  the  anisotropy  of  the  stochastic  medium,  i.e.  the  ratio  of 
the  horizontal  to  vertical  characteristic  scale,  a^/az.  The  stochastic  models  we  use  differ  in  their 
horizontal  characteristic  scales,  and  as  a  reference,  the  measurements  are  also  done  on  a  layered 
deterministic  medium.  The  measurements  are  performed  in  the  time-domain  and  take  into  account 
time-shifts  in  the  traces  caused  by  the  small  scale  velocity  fluctuations.  Under  the  single  scattering 
assumption,  the  lateral  correlation  function  of  the  wavefield  depends  only  on  the  corresponding 
lateral  correlation  function  of  the  scattering  medium.  The  stochastic  model  through  which  we 
generate  the  synthetic  seismograms  for  our  study  is  shown  in  Figure  6a.  It  is  40  km  long  and  10 
km  deep,  with  a  grid  spacing  of  20  m.  The  scatterers  are  located  in  a  5  km  thick  "target"  zone, 
extending  from  a  depth  of  5  km  to  10  km.  The  characteristic  scale  ratio  (a\/az,  with  ax  horizontal 
and  az  vertical  characteristic  scale)  of  the  scatterers  within  the  target  zone  is  varied  by  changing  ax 
to  observe  the  effect  of  heterogeneity  anisotropy  on  the  lateral  coherence  measure.  A  visco-elastic 
finite  difference  scheme  was  employed  to  generate  synthetic  seismograms.  Our  results  show  that 
changes  in  the  correlation  of  the  wavefield  does  correspond  to  changes  in  the  lateral  characteristics 
of  the  media,  though  the  differences  are  not  linear  (Figure  6b).  When  the  ax/az  ratio  of  the  velocity 
models  change  from  10  to  5  (a  factor  of  2),  the  correlation  lag,  corresponding  to  a  correlation 
coefficient  of  0.5,  changes  by  a  factor  of  1.70.  When  the  ratio  (ax/ay)  changes  from  5  to  1,  the 
correlation  lag  drops  by  a  factor  of  approximately  2.20.  This  statistical  measure  was  computed 
both  in  the  shot-gather  and  common-offset  domain  and  over  time  windows  of  500,  KXX),  and 
1500  milliseconds.  The  advantage  of  the  latter  domain  are  that  normal-moveout  (NMO)  stretch  can 
be  avoided.  Both  analysis  modes  yield  similar  results  over  the  different  time  windows,  indicating 
the  robustness  of  the  measure. 

Another  measure  we  perform  to  discern  the  horizontal  characteristic  scale  of  the  medium 
involved  a  wavenumber-time  analysis.  If  the  wavefield  is  influenced  by  the  lateral  characteristics  of 
the  medium,  then  a  lateral  wavenumber  spectra  should  exhibit  its  effect.  Comparison  of  the  power 
spectra  of  the  wavefield,  after  NMO  correction,  for  our  test  models  show  that  the  corner 
wavenumber  does  vary  in  proportion  to  the  horizontal  characteristic  scale  of  the  medium  (Figure 
6c).  The  spectra  of  the  traces  from  the  media  with  ax/ay  ratio  of  10  has  a  comer  wavenumber  of 
0.224  (1/km),  while  it  is  0.430  for  the  ax/ay=5  model  and  1.510  for  the  ax/az=l  model. 
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Ultimately  we  will  apply  these  measures  to  the  reprocessed  COCORP  and  the  PASSCAL 
datasets  for  the  Basin  and  Range  to  estimate  the  lateral  characteristic  scales  directly  from  the 
seismic  data. 

CONCLUSIONS  AND  RECOMMENDATIONS  :  The  seismic  modeling  results  show  that  a 
considerable  fraction  of  the  seismic  wavefield  observed  in  conventional  seismic  refraction  profiles 
can  be  reporduced  by  scattering  from  wavelength  scale  seismic  velocity  heterogeneity.  Moho 
reflections  are  modulated  both  by  fine-scale  fluctuations  on  the  transmission  path  and  by 
reverberation  in  the  highly  heterogeneous  near-surface  layers.  The  most  pressing  problem  to 
address  is  the  estimation  of  crustal  heterogeneity  from  seismic  data  directly.  We  have  developed 
two  measures  for  estimating  the  lateral  characteristic  scale  length  from  refraction  and  unstacked 
reflection  data.  Estimation  of  seismic  velocity  fabric  from  seismic  data  directly  will  permit 
utilization  of  the  vast  global  seismic  reflection/refraction  database  to  quantitatively  classify  the 
Earth's  crust  according  to  heterogeneity  type  and  fluctuation.  The  next  issue  to  be  resolved  is  the 
influence  of  complex  overburden  on  the  wave  field  reflected  from  deeper  layers.  Lastly  we  plan  to 
model  regional  propagation  (to  -2000  km)  in  realistic  stochastic  models  of  the  crust,  using  the 
Basin  and  Range  models  as  a  starting  point. 
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Figure  1  :  Reprcx:essed  COCORP  Nevada  Line  2,  showing  midcrustal  (4-7  s  twt)  and 
Moho  reflectivity  (9- 10s  twt).  The  upper  two  seconds  of  the  data  are  muted.  Compare  to 
the  explosion  shot  record  from  Shotpoint  4B  (Figure  4a).  Reflectivity  patterns  are  similar. 
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Figure  2  :  Stochastic  mode!  for  northern  Basin  and  Range.  There  are  no  step  changes  in 
velocity  within  the  cnist  above  ilie  Moiio.  Figure  3  :  1-D  P-wave  velocity  models  derived 
from  PASSCAL  Basin  and  Range  wide-aperture  data  by  a  variety  of  re.searchers.  Note  that 
the  mean  cinstal  velocity,  the  crustal  thickness,  and  the  velocity  gradients  are  similar  in  all 
models,  hut  the  velocity  steps  within  the  crust  are  different.  We  model  the  crust  with  a  1-D 
deterministic  velocity 
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Figure  4  ;  a)  Seismic  shot  record  4B  from  the  PASSCAL  Basin  and  Range  experiment, 
recorded  at  50m  spacing.  Note  the  bright  reflectivity  from  4-7  s  twt,  and  at  the  Moho  9-10 
s  twt.  b)  Synthetic  shot  record  for  shot  4B  from  the  velocity  model  in  Figure  2.  Note 
correspondence  between  reflective  zones  in  the  field  record. 
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Figure  5  :  a)  Composite  large  offset  shot  record  from  southern  Basin  and  Range  RISC 
experiment.  Note  prominent  crustal  events,  b)  Synthetic  shot  record  from  a  stochastic 
model  consisting  of  flat  sedimentary  basins  overlying  a  crust  intruded  by  felsic  rocks, 
similar  to  that  at  the  top  of  the  Figure  2.  A  uniform  vertical  velocity  gradient  has  been  used 
in  the  cmst  below  the  sediments. 
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